Imperial/TP/94-95/39 



Winding Number Correlation Functions and Cosmic String 

Formation 

L. M. A. Bettencourt, T. S. Evans and R. J. Rivers 
The Blackett Laboratory, Imperial College, London SW7 2BZ, U.K. 

(February 1, 2008) 



Abstract 

We develop winding number correlation functions that allow us to assess the 
role of field fluctuations on vortex formation in an Abelian gauge theory. We 
compute the behavior of these correlation functions in simple circumstances 
and show how fluctuations are important in the vicinity of the phase transi- 
tion. We further show that, in our approximation, the emerging population of 
long/infinite string is produced by the classical dynamics of the fields alone, 

being essentially unaffected by field fluctuations. 
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I. INTRODUCTION 



The production of topological defects seems to be an inescapable consequence of phase 
transitions in field theories with a topologically non-trivial vacuum manifold |IL0]. This is 
observed experimentally in numerous condensed matter systems from superfluid Helium || 
and superconductors to ordinary liquid crystals M and is expected to occur in the early 
universe for most theories of Grand Unification. Defects produced at such energy scales are 
strong candidates to generate the observed large scale structure of the universe, competing 
with the fluctuations of inflationary models to match the observations H. 

In spite of their universality, the problem of predicting the details of networks of defects 
emerging from a phase transition is still unsolved, being necessarily entangled with the 
more general question of how to determine the dynamical evolution of fields throughout the 
transition. Our qualitative picture of defect production is given by the well-known Kibble 
mechanism |]]. With the U (1) theory of a complex scalar field in mind, it essentially assumes 
that the phase transition proceeds by the formation of domains, which are characterized 
by some random value of the field phase, that later coalesce to complete the transition. 
Such domains, at the time of coalescence, can, in turn, be characterized by their average 
radial length scale £, which in most applications is taken to be the correlation length as 
computed for a massive scalar field in thermal equilibrium [|IJ]. As domains coalesce phase 
gradients must be minimized and integer windings can be formed. Because of its simplicity 
this picture has great appeal and is amenable to a straight-forward approximate numerical 
implementation by the so-called Vachaspati-Vilenkin algorithm || , which in turn constitutes 
the usual starting point for numerical network evolutions. 

As it stands the Kibble mechanism is a wholly static argument, failing to answer any 
questions regarding the role of fluctuations in the fields. If important enough, such fluctua- 
tions will necessarily alter the number and distribution of defects produced by the dynamics 
of domain coalescence and could thus distort seriously the predictions solely based upon it. 
A change in initial string density and especially in the amount of long string versus that in 
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small closed loops, for example, necessarily affects the transient regime that characterizes 
a string network prior to scaling. Moreover, it is precisely the point at which thermal fluc- 
tuations cease to be able to change string configurations appreciably that determines the 
emerging defect network. This problem has been attracting great enthusiasm leading to a 
number of different approaches [0-0. 

Our aim in this paper is to calculate the effect of fluctuations on winding number, 
requiring us to develop quantities and criteria for describing fluctuations about general 
field configuration solutions to the phase transition dynamics. We focus our analysis on 
cosmic strings in the Abelian Higgs model and consider correlation functions of the winding 
number (defined as the number of quantized units of magnetic flux) in several circumstances. 
Even though the present study is motivated by cosmological questions the analogy with 
superconducting systems is almost complete, because of the role of the gauge fields in both 
models. As such our results should apply equally well there, where contact with experiment 
can prove a good test on the ideas developed herein. In principle the interpolation between 
both pictures can be obtained essentially by replacing cosmic strings with vortex lines and 
discarding some specific aspects of relativistic self-energies in section 3. In practice we have 
yet to so so. 

The outline of the paper is as follows. We start, in section 2, by describing the qualitative 
aspects of the classical dynamics of the fields during the phase transition and determine under 
which conditions the measurement of quantized units of magnetic flux becomes a reliable 
estimate of the number of strings present. We then proceed, in section 3, to compute 
the fluctuations in winding number in the simplest scenario, namely that of a constant 
background scalar field with an equilibrium distribution and under the assumption that the 
gauge field, too, displays a thermal spectrum. In section 4 we show that the computations 
of section 3 can easily be extended to the calculation of correlation functions of the winding 
number threading through two different contours with relative specific geometries. We use 
these correlation functions to answer questions concerning the conditional probability of 
finding string crossing a given surface in space, subject to the condition it has crossed 
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another, with a given separation and orientation relative to the former. In section 5, we 
return to the computation of winding number fluctuations, by including inhomogeneities in 
the scalar field. We show how a systematic expansion around a homogeneous background 
can be performed, generating as a zeroth order term the result of section 3, and compute 
the first order correction resulting from the presence of string. By comparing their relative 
magnitude we can formulate a well defined criteria for the breakdown of such an expansion. 
In section 6 we present our conclusions and compare our results with those of other recent 
approaches to the same problem. 



II. PHASE TRANSITION FIELD DYNAMICS AND DEFECTS 



We begin with the semiclassical analysis of the simplest theory that admits gauged 
strings, the Abelian Higgs model. We consider it in the symmetry broken phase so that the 
decomposition of the complex scalar field into modulus ip and phase a is possible as 



</>{x) =p{x)e~ ia{x) . 



(1) 



for non-zero (fi(x). Thus the classical action, including the gauge fixing term, can be written 

as 



-4V+2 



+ —d^A^d u A u . 



(2) 



Our objective is to compute string formation in gauge theories, as a dynamical process 
occurring as a by-product of the symmetry breaking in the theory, by measuring units of 
the quantized magnetic flux crossing a surface bounded by the contour 7, given by 

(N 7 ) = ±£dT(A i ). (3) 

The idea of using the magnetic flux as an indicator of the net winding number springs 
from the vacuum field configurations in the broken phase, namely, 



(<p) = V (4) 
(A, ~ d -f) = W = 0. (5) 

This guarantees that the gauge field indeed traces the gradient in phase of the scalar field. 
In a situation where the scalar modulus is constant, however, this becomes simply a gauge 
transformation, corresponding to the residual freedom in defining in the Lorentz gauge 
(d^A 1 * = 0). However, it is the case of a non-constant scalar modulus that is relevant for 
our purposes, corresponding to the presence of strings. In this case the phase can change 
around a closed contour in space by a non-zero multiple of 2it, implying the existence of a 
singularity in its gradient as the contour is shrunk to a point (at which there is a necessary 
restoration of symmetry). The gauge field at such a point, the core of the string, exhibits 
therefore its expectation value in the symmetric phase, i.e., (A^) = 0. The approximate 
solution of the static Euler-Lagrange equation interpolating between these two regimes is 



well-known to be 10 



Ti 

A e = n(pKi(e(pr), (6) 

er 

for a homogeneous scalar modulus (p. At large distances from the the core of the string ( 
r < ), the modified Bessel function K\ becomes exponentially decreasing and the vacuum 
of the spontaneously broken theory is rapidly approached. In this limit, applying the form 
(6) in (3) gives precisely n units of flux with exponentially small corrections. 

As the phase transition proceeds and the fields become smooth we therefore expect A^ 
to trace the overall differences in phase left over from the dynamics of the scalar field with 
increasing accuracy. To see in more qualitative detail how this happens we need to consider 
the actual Euler-Lagrange equations of the theory 

A 



d^ip = if 



(eA^ - d,a){eA» - d»a) - '-{ V 2 - V 2 ) 



(7) 



3" ( + -d u A^ ) = J„, (8) 



a 



and 
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<9% = 0, (9) 

where 

J M = ey? 2 (e,4 M -c». (10) 

The phase transition is triggered by an instability in the long wave-length modes of the 
scalar fieldP], whose amplitudes then start growing exponentially towards r\. This is expected 
to happen independently in regions of space-time which are causally disconnected for the 
time scales involved. As a result a domain structure is expected to form, leading to an 
inhomogeneous scalar modulus on a scale of several domains. These inhomogeneous field 
configurations do not minimize energy and will therefore evolve in time so as to coarsen and 
smooth out 0. 

What is the behavior of the phase and vector potential during this dynamical process? 
The vector potential is characterized initially, in the symmetric phase, by (A^) = 0. During 
the transition A^ is sensitive to the change in expectation values of cp and a in quite a 
different manner. Whereas the change in the expectation value of tp gives a mass to A^, 
without changing its expectation value, the presence of a gradient in phase does act as a 
source for the latter. This can be seen by writing (8) as 

[d^A„ - (1 - -)d u d»A,] - e 2 <p 2 A u = - ecp 2 d v a = J'„ (11) 

(X 

where J^ xt is any external source, imposed on the system. 

The general (A^) can, in principle, be computed from the generating functional for the 
theory which in turn follows from expanding the action (2) around a field configuration that 

1 The details of the triggering mechanism depend crucially on the hierarchy of the couplings in 
the theory determining thereby the order of the transition, as seen from thermodynamic effective 
potential constructions. 

2 This is explicitly observed in condensed matter systems, such as specific types of liquid crystals, 
which display a large relaxation time Q. 
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satisfies the Euler-Lagrange equations with boundary conditions characterizing the fields 
before and after the phase transition, as discussed above. 
By introducing (11) in the action, we obtain 

1 t 

S[ip, a, Afj\ = -- J d^x (e^d^a - J^ xt ) + terms in ip and a 

= ^Jd 4 x J'^ + terms in ip and a . (12) 

This can in turn be re-written exclusively in terms of the currents, as usual, since 

A lt (x) = Jd 4 yG^(x,y,cp)j' l/ (y), (13) 

where 

d»d„ - (1 - -)d u d» - e 2 <p 2 (x)] G, u (x,y,<p) = -8 v (x - y). (14) 

In the absence of external sources A^ then becomes 

A ix( x ) = J dA y G fM/ (x,y,(p)eip(y)d' / a(y). (15) 

As the phase transition approaches completion the scalar field should become smoother, 
both in modulus and phase, leading to small derivative terms relative to the mass scale. 
This allows us to perform a small derivative expansion of the gauge propagator, leading to 
the result 

. . . / (p(x) 2 \ dna(x) . . . . 

AJx) ~ — — — h derivative terms m ip and (16) 

\ i] 2 J e 

In the case of a homogeneous and static scalar modulus, this reduces to the vacuum in the 
spontaneously broken phase (5). Unlike its higher derivative corrections, the first term in 
(16) is gauge independent. 

Having understood qualitatively from this semi-classical analysis how the gauge field 
reacts to a difference in phase, it is important to discuss how a gradient in phase itself can 
occur and how it will evolve. Departing from the symmetric phase, where (A^) = 0, eq.(9) 
can be written as 



d lx d fX a + (<9 M ln</? 2 ) - 



(17) 



Eq (17) is homogeneous in d^a showing that a gradient in phase must exist initially for the 
phase to change in space-time. Simultaneously, given that = 0, the energy of a given 
configuration is minimized for a constant value of a. We therefore expect the field phase to be 
approximately constant inside each domain formed by the dynamics of the scalar modulus. 
Differences in phase can then only occur when different domains come into causal contact, 
and then should evolve so as to be minimized. This picture corresponds to the well known 
Kibble mechanism, in a gauge theoryfj. In particular, according to (17) and under plane 
wave ansatze for the fields, the small frequency modes of a will be resonant with those of ip. 
As a result, given an exponential growth in the amplitudes of ip a corresponding exponential 
decrease in a will occur. As domains coalesce and d^ip tends to zero the remaining gradients 
in phase are frozen in and the gauge field catches up with them, in the sense of (16). 

Having this qualitative picture in mind for the description of the field dynamics during 
the phase transition we proceed, in the next sections, to compute the effect of fluctuations 
on the total winding number using the quantized magnetic flux (3) as a probe. These are 
given by 



In most of what follows we will be concerned solely with thermal fluctuations. Quantum 
fluctuations can be readily included but are known to be subleading. Their effect is to 
change the value of the vacuum expectation value of the yj-field as well as the couplings 
in the theory. We will in general adopt the position of treating <p as providing a classical 
temperature dependent background field (<p) in which the fields can be quantized. In general 

3 We see that the geodesic rule comes about naturally from dynamical considerations. Ambiguities 
PI result from considering the static gauge theory in its spontaneously broken phase, where the 
phase decouples and loses physical meaning. 




(18) 
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this background will include strings. We then need to know the form of the two point function 
for the gauge fields, which in the absence of external sources, is given by 



where G^ v is given by (14), and J M (x) = eip(x)d^a(x). If we assume that the contours do 
not run over any strings^ the dominant term will be 



We stress that we are not just quantizing A^. Scalar-field radiative corrections about (<p) 
will be taken into account in the calculation of G^. The remaining uncertainty in defining 
the gauge propagator explicitly results from the unknown form of <p{x). Having (16) in mind 
we will treat the situation of a quasi-homogeneous background, expanding the propagator 
in the perturbations. The time dependence of ip remains a problem, however. We deal with 
it approximately by considering the fields in equilibrium and changing the corresponding 
temperature. In the next section we treat the simplest case, i.e., that of a homogeneous 
background where the fields display equilibrium distributions. 

III. FLUCTUATIONS IN WINDING NUMBER IN THERMAL EQUILIBRIUM 

In this section we specialize the previous discussion to the case of a homogeneous scalar 
modulus background and assume the fields to have equilibrium distributions. The first of 
these conditions will be lifted later, in section 5. Dropping the condition of equilibrium, 
however, would require the knowledge of the time dependence of the fields, in closed form, 
which is unknown in the necessary detailf]. 

4 This seems reasonable if the string density is not too high. See section 5. 

5 More is known for the theory of a simple scalar field numerically, e.g. |12], but this proves of 
little use in our present context. 



{TA^{x)A v {y)) = G^x, y, <p) + {J"{z)G m {x, z, v)r(z')G av (z\ y, <p)) 



(19) 




(20) 
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The degree of approximation actually involved in adopting a thermal distribution for the 
fields will necessarily depend on the dynamics of the phase transition, which in turn will 
be determined by the hierarchy of couplings in the theory. Firstly, for a theory with strong 
gauge couplings (Type I) one expects the system to depart significantly from equilibrium 
(by bubble nucleation) and our approximation will probably grossly misrepresent the actual 
fluctuations in the system. In the opposite limit, however, (strong Type II) it should hold 
as a good approximation more or less throughout the whole dynamics. The latter is the 
case of most high temperature superconducting systems. Secondly, in any case, the late 
stages of the transition should be characterized by the approach to a thermal distribution 
and consequently the estimates below should describe the actual fluctuations in the system 
increasingly more accurately. This expectation is also consistent with the indications of 
when the magnetic flux becomes a good indicator of the number of strings crossing a given 
surface, as seen in the last section. 

A strictly homogeneous scalar field is of course a necessary and sufficient condition for 
the existence of no strings. It then becomes necessary to be able to estimate how the values 
of the fluctuations computed below relate to those in more realistic inhomogeneous back- 
grounds. This issue will be addressed in section 5, where we will show that under well-defined 
circumstances the fluctuations computed here are the first term in a systematic expansion 
about a homogeneous background and constitute, in general, the dominant contribution. 

In order to compute the thermal fluctuations in winding number, given by (18), we need 
to know the form of the time-ordered thermal correlation function of the A^. There are 
only two symmetric tensors in the spatial indices which can be made in the covariant and 

gauges, namely 5 lj , k l k j . The Ward or BRS identities and the Lorentz structure at finite 
temperature then tell us that the most general form in covariant and gauges is 



iG ij (t -t',x- x') = (TA i (t, x)A j (t', x')) 



(21) 




(22) 
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-G ij (k) 



1 

k 2 — e 2 r] 2 — Ut 
( 1 

+ 



(6 



ij 



k l W . 



+ f(k,M,() 



k l ti 



(23) 



x fc 2 — e 2 rf — Hl ' ' J k 2 

where the energy fc = 2irin/ f3. The / function varies with the gauge chosen while the other 
il-r and Hz, terms correspond to the physical modes of the photon, two transverse (magnetic) 
and a single longitudinal (electric) mode. 

From (18) and (21-23), taking the integrals over circular paths of radius L, we see that 
we require two sorts of integration 



dxie 



-ik.x 



do£j& k 



(24) 



Since the latter is zero we only pick up the term in the propagator with the Sij. This is not 
surprising as we are really looking at the magnetic field correlations and the term with 8^ is 
associated with magnetic fields, its self-energy, ITy contains information about the magnetic 
screening. Thus we have that 



E 



d A k 



—% 



d>JCj 6 



-ik.x 



n{uj) 



(25) 



4tt 2 8tt 3 /3 n f^ 00 J k 2 - e V - n T 

Note that this form is independent of the gauge chosen and the only source of gauge depen- 
dence is in the expression used for Hy. Proceeding with (25), we use the leading term in the 
high temperature expansion of the one-loop self-energy. In this we are working to leading 



order in the resummation scheme of Braaten and Pisarski |T3| which is necessary at high 
temperatures. Performing the energy sum leads to two terms. There is a pole contribution 
which corresponds to propagation of a 'plasmon'. The pole is found at ko = ±a>(fc) where 



= u{kf - e 2 r](T) 2 - IL T (u(k), k) 



(26) 



n T (k ,k) 



3M 2 fh 



k 2 [ k 2 2k 



ko — k 
k + k 



where we have worked with respect to the vacuum 



0. 



(27) 
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The usual 0(e 2 ) result is rj(T) 2 = rf (l — which vanishes at the critical temperature 
T c . For abelian theories, the value of M p is 

Ml = 6 -j-, (28) 

for a U(l) theory with a single charged scalar or fermion. While the position of the pole for 
general three-momentum is given by a non-trivial dispersion relation, to within about 10% 
this is approximately expressible as k 2 , = k 2 + m 2 where 

m 2 = e 2 r](T) 2 + M 2 P . (29) 

The difference in the final results, between choosing this approximate dispersion relation or 
the true one, has been studied numerically in some of the cases below and is found to be 
negligible (less than 1% ). 

The real complication is encountered with the cut term, running across the zero energy 
point, seen in the imaginary part of Ht for |fco| < k. This corresponds to the physical process 
of Landau damping allowed in a heat bath. This must be treated numerically as in this case 



we are not able to use any sum rules of ||14|| , which can be seen in the logarithmic dependence 
on temperature. Equivalently, trying to look at possible 0(T) terms by including only the 
n = term in the energy sum (e.g. see appendix of [fL5|| ) leads to divergent integrals and to 
the breakdown of the calculation. 
Doing the energy sum gives 



(at 2 ) = -i- 1 d 3 k 



dx^e lk ' x 

L 



2 ( n(uS) 



+ (cut terms) 



\ uj J 
+ (T=0 terms). (30) 

On computing the integral in (25) we find 

(N 2 ) = -^L 2 J k\\dk\\dk z Jl(k\\L) (^^- + (cut terms)^) + (T=0 terms), (31) 

where k\\ is the momentum component on the plane defined by the contours and J\ is the 
usual Bessel function of order 1. This last integral cannot be fully computed analytically. 
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For large enough loops (k\\L < 1) we can take the large argument form of the Bessel function. 
Looking only at the pole contribution and using the good approximation that oj 2 = k 2 + Mp 
leads to 

(AO-^^lnf-V (32) 



47T 2 3 vm, 

This estimate is indeed confirmed by computing the integral (31), for the pole term, numer- 
ically. Figures 1 and 2 show the result of a numerical integration looking at just the pole 
term with the approximate dispersion relation oj 2 = k 2 + m 2 on the dependence of (N 2 ) 
on LT and respectively. The ~ coefficient results from the best-fit, given the functional 
form (32). Using the exact dispersion relation for the pole contribution makes only a very 
slight difference to the results. The cut term shows a similar sensitivity on T](T) which acts 
as an infra-red cutoff. For Lrj(T) < 1 the cut contribution to N 2 rises as LlnL. However, 
for Lr](T) ~ 1 and larger, a similar qualitative behavior to (32) is found, with magnetic 
screening mass [] m mag = erj(T) replacing the plasmon mass in the formula above to give 

<^ 2 > « ^ LThl (w))- (33) 

with a coefficient of proportionality comparable (typically, a few times larger) to that of the 
pole in (32). Thus we see that the contributions from the pole and cut terms essentially 
generate the same form for the winding number fluctuations away from the critical temper- 
ature, when T)(T) <C Mp. As such the influence of the cut term, in these circumstances, is 
to change the overall coefficient in the functional form (32). In most of what follows we will 
therefore restrict ourselves to computing pole contributions. 

Qualitatively the results above can be thought of as follows. The dependence of (N 2 ) on 
L can be understood in terms of a random walk in the field fluctuations along the contour 
perimeter. This is commensurate with the Kibble mechanism of domain formation, with 
random phases correlated on a scale £ = l/e 2 T (up to logarithms). 



6 The magnetic screening mass for abelian theories only has a contribution from the scalar field 
shift. 



13 



The quasi-linear dependence on the temperature ensures that, as the system cools, fluc- 
tuations in winding number become less and less relevant. The logarithmic dependence on 
the field thermal mass, in turn, guarantees that as one approaches the critical temperature 
and the mean field mass r\ (T) vanishes, the magnitude of winding number fluctuations di- 
verges. This is certainly the instance when thermal fluctuations can significantly change 
any underlying expected winding number. The actual divergence of (N 2 ) at the critical 
temperature for the pole term, however, is an artefact of our best fit since, as can be seen in 
fig. 2, as the mass vanishes (for Mp ~ 0) the fluctuations in winding number remain finite 
even though maximal. This is not the case for the cut contribution, as far as we could find 
numerically where the logarithmic behavior persists, with extraordinary accuracy, down to 
m = IfT 4 

'"rnag ^ u 

By themselves, however, the magnitude of fluctuations tell us little about how string 
configurations can be changed. This then takes us to the crucial issue of comparing the 
magnitude of fluctuations to the underlying expected winding number. The magnitude of 
relative fluctuations is given by 

(34) 

(TV) 2 (iV> 2 

where (A} is the total string number, without distinguishing strings from antistrings. In 
principle (A^ 2 ) and (N) 2 must be computed self-consistently given a specific field background. 
To estimate the behavior of (34) in our specific setting, however, we can take our values for 
(A^ 2 ) and borrow the value for (A^) 2 , obtained using the Vachaspati-Vilenkin algorithm, for 
the number of strings crossing a disc of area nL 2 . As required, this algorithm does not 
distinguish between strings and anti-strings. Using this value in (34) should clearly yield an 
underestimate for the value of the relative fluctuations. The expected value of the winding 
number thus computed crossing a disc of radius L is f| 

<*> = ~ (35) 
where £ is the average domain radius at coalescence. This can be written as an inverse of a 
mass scale £ = — , whence 
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(AiV) 2 8e 2 T 1 




(36) 



(TV) 2 3tt 2 m^L 



The identification of this scale with the mean domain radius at coalescence springs directly 
from its numerical implementation. Its computation from the field theory, however, implies 



proposal by Kibble [|TJ was to identify this scale with the inverse Ginzburg temperature 
which is a well defined quantity in a field theory undergoing a second order phase transition. 
More recently, Kibble and Vilenkin argued that a more suitable scale could be found by 
invoking the scaling of the string network emerging from the phase transition. They discuss 
carefully the different length scales in the problem and essentially conclude that individual 
strings then can be identified once this typical string separation scale becomes larger than 
the Ginzburg length. Before then fluctuations in the fields are argued to be very large 
and coherent long-lived field configurations cannot persist. Looking at (36) we can indeed 
observe that fluctuations become large for temperatures larger than and length scales 
smaller than its inverse. At the critical temperature the relative fluctuations should diverge 
due to the cut contribution discussed above. 

Our present analysis thus lends support to several of their qualitative arguments and 
clarifies the emergence of some scales from field theory propagators. 

IV. STRING SPATIAL DISTRIBUTION AND OTHER CORRELATION 

FUNCTIONS 

One interesting question to ask about string formation that is of the upmost importance 
for the subsequent evolution of a string network is what is the fraction of string in loops as 
opposed to long string. We will see in this section that winding number correlation functions 
can give us some insights into this problem. 

Since we expect to be dealing with a complicated background in which many string 
configurations will thread through our contour these questions can only be answered statis- 
tically. Intuitively then, the question of how string is distributed in space assumes the form 



a knowledge of the domain dynamics which is, in general, very poor |16|]. The original 
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of a conditional probability question, i.e, knowing the amount of string that crosses a given 
surface what is the amount of string that crosses a second surface with a given orientation 
relative to the former. Answers to these questions are naturally provided in terms of cor- 
relation functions and the natural way ahead, in the context of our present discussion, is 
to take our two point function for the winding number (18) through two different contours. 
In so doing we expect that the correlations between net winding number threading through 
the two loops in a given geometrical configuration can also give us information about how 
string is distributed relative to that geometry. The general case is too difficult, but there 
are two natural extensions of (25) that are solvable, namely that for which the two circles 
remain concentric but now have different radii / and L > I and that of two co-axial circles 
of equal radius L separated by a distance h. These are shown in Figs. 3 and 4, respectively. 
Intuitively we would expect the first of these correlation functions to give us information 
about the population of string with curvature radius greater than that of the longest con- 
tour whereas the second should provide a measure of straightness on the scale of separation 
between the contours. 

We start with the former. In what follows we compute the contributions to the correlation 
functions due to the pole term only. As discussed above the cut contribution should not 
change any qualitative aspect of our results but merely the overall coefficient by a number 
of order 1. As before in order to compute the the two-point function for the winding number 
we can show that only the integrals of the type I\ survive, i.e., we have to compute 

h = J dlidlje^-^Sij = IL J" ded6'e ik ^ Lcosd '- lcos V cos(0' - 6). (37) 

This can be shown to yield 

h = 4tt 2 /LJ 1 (A;||L)J 1 (A;||/). (38) 

We then see that in this case the correlation function is no longer, necessarily, positive 
definite but rather depends on the interplay between the radius of the two contours. The 
expression for the correlation function then becomes 
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( NN ) = * U I ^dfc. jlMMM fl^), (39) 



4tt 2 ' 

Figure 5 shows the result of the numerical integration of (39), having fixed the radius of one 
of the contours to be I — ^, ^ and ^ and m = 0.2T. It shows that the winding number 
is most correlated when the radius of the two contours is comparable. Moreover, when the 
difference between the radii of the two contours becomes appreciable the correlation function 
exhibits the approximate behavior 

2 -1.25m(L-l) 

(N 7 Ny) ~ — LT0.05— (40) 

Thus, for radial separations larger than essentially a few units of the inverse mass, the 
winding number through the two contours becomes essentially uncorrelated showing that it 
should correspond to different strings. String, if produced by a thermal fluctuation, would 
therefore be expected to have a curvature radius smaller than that scale, corresponding 
necessarily to loops of that form. 

Another possible generalization of (25) corresponds to translating the two contours, in 
the direction perpendicular to the plane they define, relative to each other so that they now 
define the two basis of a cylinder of height h, see fig.4. In this case the form of the original 
integral still holds but the Fourier transform corresponding to the V-direction ceases to be 
trivial. We then obtain 

(N^Ny) = ^-L 2 J k\\dk\\J?(k\\L) J dk z e ik * h -^, (41) 

where h is the separation between the two contours along the vertical axis. This gives 

(ALiVy) = 2 — L 2 / h\dh\J 2 Ah\L) dk z cos(k z h)^^. (42) 
An Jo Jo to 

Figure 6 shows the result of these integrations for several values of L and m = 0.2T. Again 
an approximate analytical behavior can be found for (42). The best fit on the mass and h 
dependences yields 

e 2 

(NyNy) ~ — 0.6LTe- L25mA (43) 
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Both correlation functions (39) and (42), above exhibit exponential fall-offs with contour 
separation, determined by a correlation length multiple of the inverse gauge field mass. This 
behavior strictly holds for m(L — I), mh))l. Together with the result of the previous section 
about the magnitude of fluctuations, these allow us to characterize the temperature for 
which winding number fluctuations are relevant on a given length scale that can in turn be 
compared to the string's proper length. In is clear that, as the temperature drops, only the 
population of string loops smaller than £ = O(^fy) will be affected. Thus any long/infinite 
string emerging from the phase transition will be very little changed by fluctuations. This 
is an important result as it states that the formation of long string is essentially a classical 
process occurring at the level of the Euler-Lagrange equations for the theory. 

V. FLUCTUATIONS IN WINDING NUMBER IN THE PRESENCE OF STRINGS 

In this section we investigate how a non-homogeneous scalar modulus changes the fluc- 
tuations in winding number, as computed in section 3. This is an obvious requirement 
of any description of fluctuations taking place during a phase transition. Inhomogeneities 
occur naturally in the dynamics of domain growth and coalescence, corresponding to such 
interfaces as well as to topological defects}]. As the phase transition approaches completion 
we expect the domain structure to coarsen and smooth away and the only inhomogeneities 
left behind should correspond to topological defects. At all stages, thermal fluctuations in 
winding number will occur and a freeze-out in defect number will correspond to the instant 
in time when their role is rendered negligible, in some well-specified way. 

In what follows we will attempt to give a general treatment of fluctuations in the pres- 
ence of a weakly inhomogeneous scalar modulus by systematically expanding the gauge field 
propagator in a small perturbation around a homogeneous scalar field. We will later special- 
ize to the case where strings constitute such inhomogeneities. As before we will only keep 



7 In the early stages of this dynamics such a distinction probably fails to apply. 



18 



the pole term in the gauge field propagator. 

Let us consider the case of a localized deviation from the otherwise homogeneous back- 
ground so that we can write the scalar modulus as 

(p(x) = 7] - f(x) (44) 

where 77 is the background value of the scalar field modulus and f(x) the small perturbation. 
Then we can expand the gauge propagator around the homogeneous background field. We 
restrict ourselves to the relevant part of the transverse propagator, for the sake of clarity. 
Formally, we obtain 



G^(x,y,M) = (xl-j- rTT72? \ \v) ~ ( x \ n2 ^ 2 

O z + m z + M z (z) D z + m z 



l-\z)M 2 (z)(z\ 1 2 + 
u z + m z 



+ \\z)M 2 {z)(z\ - 1 - \z')M 2 (z')(z'\ - 1 - 
2 1 ' v M '□ 2 + m 2 ' Q 2 + m 2 

+ higher order terms] (45) 

where M 2 {x) = e 2 f(x) [2rj — f(x)]. The first term in the expansion is just the gauge field 
propagator in the presence of a homogeneous, static scalar field in the broken phase. When 
used to compute winding number correlation functions it will generate the results of the 
previous two sections, when taken in thermal equilibrium. The first correction to that can 
be written as 

G$(x, y) = 8^ J d 4 pd 4 kd 4 zd 4 z'G{p)e ip(z -^M 2 {z, z')G{k)e lk(y - z '\ (46) 
where M 2 (z, z') = M 2 (z')5\z - z'). This leads to 

G$(x, y) = 5, u J d 4 pd A kG(p)e- ipx M 2 (k - p)G(k)e lk y, (47) 

where 

M 2 (k-p) = J d 4 ze- iz ^ k - p) M 2 (z), (48) 

is the usual Fourier transform of M 2 (z). We see that the presence of inhomogeneities can 
be accounted for as the occurrence of sources, in the usual way. Their effect is naturally 

19 



to change the 4-momentum of photons from an in-going state to an out-going one. The 
order in the expansion (45) then corresponds to the number of such momentum transfers. 
To actually compute corrections to our previous results we have to specify the functional 
form of M 2 (x). As an illustration of the effect of inhomogeneities we take the simplest case 
of associating M 2 with one straight string at rest, placed in the center of the loop. That is, 

M 2 (x) = e 2 f(x) [2r] - f(x)} ~ 2e 2 T]f(x), (49) 

where we specialize the form of f(x) to correspond the the field profiles of a string at 
distances larger than its width [lTj 



/(r) = k s K (m s r) ~ J—I_ e -"^, (50) 

where Kq is the modified Bessel Function of order zero and ks is a constant that for critically 
coupled strings can be shown to be ks = \n\rj, with n as the string's winding number. Here we 
assume strings to be the usual static solutions of the effective 3-dimensional theory, obtained 
after the usual frequency sums at finite temperature. As such all mass scales exhibit their 
usual temperature dependence. The exponential regime corresponds to the asymptotic form 
of the Bessel function for large arguments. In particular we see that, in this regime, the 
corrections to the homogeneous background propagator will be exponentially suppressed by 
the distance to the string. 

It is straightforward to generalize the situation for one string to an arbitrary number of 
strings, or other well localized sources of inhomogeneity with a definite profile. We simply 
write the scalar modulus as the superposition of all sources of inhomogeneity as 

N 

if = r] 1 ~ N Ilf =1 (N — fi) ~ 1] — fi + higher order terms in f (51) 

i=i 

Let us then proceed to compute the change in fluctuations brought about by a single string, 
to first order. The form (50), corresponding to a single static cylindrically-symmetric string 
placed in the center of our circular contours, actually allows us to compute its Fourier 
transform analytically; namely 
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M 2 {k-p) ~ J d 4 xe' l(p - q)x 2m 2 e K {m s r) 



= ^m f*-*M k '- p '\ (52) 
(k — p)j + mj 

which is just the potential, in momentum space, for a massive cylindrically symmetric field, 
as it should be. Using this form in (47) yields the correction to the winding number corre- 
lation functions 

Ji(p||L)Ji(fc||L) / v\ + k\ + m 2 



7 8^ e} J p W 11 Wq w k ^(pj + jfcj| + m2)2-4pffc 
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X 



(n(u> p ) + n(iwfc)) (n(w k ) - n(w p )) 



(53) 

_w p + w k x w k -Wp x 
where we have assumed that both the in-going and out-going momenta display thermal 
distributions. This follows automatically since a static source causes no energy shift. 

The integral (53) has to be computed numerically. Its dependence on L, together with 
that of the zeroth order result of section 2 and their sum are depicted in fig. 7, for ^ = 
0.5, ^ = 0.6. It shows that the effect of this correction is only to change winding number 
fluctuations significantly for small contour radius, smaller or equal to the string's width. 
This can also be seen by fixing L and m e and varying m s . Indeed, the magnitude of the 
correction term becomes larger the smaller scalar masses we take, i.e., with greater string 
widths, signaling the breakdown of the weakly inhomogeneous scalar field expansion for 
the gauge propagator. In this regime the sum of the two terms yields a negative quantity 
showing the clear breakdown of the expansion (45) since the full quantity is positive definite. 
For large contour radius the effect of the correction becomes increasingly more negligible 
and the negative shift probably becomes a genuine feature. 

Another relevant issue is to know how this correction term behaves close to the critical 
temperature. This can be probed by assuming that both masses vanish approximately 
as we approach the critical temperature. In practice we know this to be untrue for the 
gauge field due to higher order effects, as discussed in section 3. We can distinguish two 
cases, namely, when the Yukawa gauge coupling is larger or smaller than the scalar self- 
coupling. This differentiates between a type I and a type II theory, respectively. Fig 8, 
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shows the dependence of the zeroth order term on the gauge mass as well as that of (54) 
for m s = 1.5m e , m e , 0.5m e , for fixed L = 5/T. We see that the magnitude of the correction 
term grows with the gauge coupling, while still remaining smaller than the zeroth order term. 
This however does not hold for strong type I theories for which we see that larger values of 
L will be required for the series expansion to hold. One important feature of the correction 
term, that holds regardless of the hierarchy of the couplings in the theory, is that it vanishes, 
at critical temperature, with the scalar mass. This should constitute a genuine behavior of 
all correction terms since they are proportional to increasing powers of M(x) . Thus as 
the latter goes to zero so does M(x) and we recover the usual electromagnetic propagator. 
This allows to conclude that the weakly inhomogeneous background expansion can always 
be trusted in a sufficiently small neighborhood of the critical temperature and that the effect 
of the presence of strings will manifest itself mostly at intermediate and low temperatures. 
These conclusions would radically change if we do not adopt a thermal equilibrium mass 
scale (vanishing at the critical temperature) in our defect configurations. In assuming this 
we hope to mimic an average defect configuration in the presence of a thermal plasma and 
thereby render the whole picture self-consistent. 



VI. CONCLUSIONS 

In this paper we have developed techniques to allow us to study the role of fluctuations in 
cosmic string formation in an Abelian gauge theory. Their application to the simplest case of 
quasi-homogeneous scalar field backgrounds leads us to conclude that thermal fluctuations 
in winding number are strongest the closer the system is to the critical temperature and, 
further, that the relative winding number fluctuations diverge at that temperature signaling 
the presence of no stable defects, in these circumstances. This effect persists even in the 
presence of inhomogeneities in the fields, such as strings themselves, if thermal equilibrium 
is used throughout. As the temperature drops the role of thermal fluctuations in changing 
winding number becomes less and less important. The energy scale that characterizes an 
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approximate freeze-out in defect numbers is most naturally the inverse mean domain size 
at the time of coalescence, which in turn is associated to the average string interseparation. 
From the point of view of the field theory its determination necessarily implies a more 
detailed knowledge of the field dynamics during the transition than what is currently known. 

On developing generalizations of our initial winding number correlation function we also 
conclude that string configurations can only be changed by thermal fluctuations on scales 
of the order of the thermal correlation length of the fields. As such the population of small 
string loops can be drastically modified but not that of long and infinite string, which is 
the most relevant input for network numerical evolutions. A change in loop population is 
only likely, for the present cosmic string network implementations, to change the approach 
to the scaling regime and will therefore bring no change to the details of late time structure 
seeding in the canonical scenario. 

An analogous winding number correlation function over the same contour can be defined 
using the phase of the Higgs field directly 0. Using the gauge sector instead has great 
advantages in computational simplicity and reliability as the former is divergent as the 
symmetric phase is approached from below. This leads to infinities and computability can 
only be maintained by on adopting severe approximations. Both methods will necessarily 
measure fluctuations in the gauge and scalar fields, respectively, that do not necessarily 
correspond to strings. This is because, even though one of the fields can have the right 
configuration to generate a string, it will take the two at the same point in space to produce 
a Nielsen-Oleson vortex. One is then overcounting strings as we will be computing the 
number of strings present at a given instant together with that of fluctuations that can 
potentially produce a string or decay away 0. 



Other criteria exist to count the number of strings and their fluctuations, namely those inspired 
in counting the density of zeros of the scalar field IfTql . However, the same methods have yet to be 
applied to gauge theories. 
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Our methods should have other applications beyond the GUT transitions. The system 
considered above has obvious analogies with a superconductor where the Ginzburg-Landau 
free energy is given precisely by a non-relativistic analogue of the Abelian Higgs model. Our 
analysis for the winding number correlation functions generalizes to both models provided 
the mean field masses and self-energies get properly translated. This has yet to be done, but 
would possibly constitute the natural test on the approach developed above as comparison 
with experiment seems a feasible task. The next closest example of string configurations 
in a specific model concerns those in the Electroweak Standard Model. Their study in the 
simplest approach consists in treating the model as an effective Abelian theory of a complex 
scalar field and the Z gauge boson. Even though the supplementary field excitations are 



known render them unstable |T9[ our winding number correlation functions should find 
applicability here, too. 

The eventual applicability of this approach in genuinely out-of-equilibrium scenarios 
hinges upon a better knowledge of the corresponding fields and their two-point correlation 
functions. A considerable amount of effort is being devoted to this problem JT2],^U| and it is 
conceivable that the full problem will be tractable in the near future. 
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FIGURE CAPTIONS 



Figure 1 : The dependence of the winding number fluctuations on the contour radius 
L from the pole contribution, for m = 0.2T, 0.5T and 0.8T. 

Figure 2 : The dependence of the winding number fluctuations on the gauge field mass 
from the pole contribution, for L = T, 5T and 15T. 

Figure 3 : A Schematic view of the two concentric contours and of the strings that 
contribute to the winding number correlation function between them. 

Figure 4 : The same as in Figure 3 for the two coaxial circular contours. 

Figure 5 : The dependence of the winding number correlation function between the two 
concentric contours on one of the contour's radius L, for m = 0.2T and / = T, 5T and 10T. 

Figure 6 : The different curves exhibit the dependence of the winding number correlation 
function between the two circular coaxial contours on their separation ft,, for m = 0.2T and 
L = 5T, 10T and 15T. 

Figure 7 The magnitude of the zeroth order term (dotted line), the absolute value of 
its first correction and their sum as a function of contour radius and for mg = OAT and 
m = 0.2T. 

Figure 8 : The dependence of the zeroth order term in winding number fluctuations 
(dotted line) and its first order correction on the gauge field mass M, for L = 5/T and in 
the cases mj = 0.7m, m and 1.2m. The mass scale M is in units of temperature. 
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